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ABSTRACT 

We investigate the origin and nature of the profiles of water and ammonia 
observed toward the L1157-B1 clump as part of the HIFI CHESS survey (Cec- 
careUi et al. 2010) using a new code couphng a gas-grain chemical model with 
a parametric shock model. First results from the unbiased survey (Lefioch et 
al. 2010, Codella et al. 2010) reveal different molecular components at different 
excitation conditions coexisting in the Bl bow shock structure, with NH3, H2CO 
and CH3OH emitting only at relatively low outflow velocities whereas H2O shows 
bright emission at high velocities. Our model suggests that these differences are 
purely chemical and can be explained by the presence of a C-type shock whose 
maximum temperature must be close to 4000 K along the Bl clump. 

Subject headings: ISM: individual objects: L1157 — ISM: molecules — stars: 
formation 
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1. Introduction 

During the earliest stage of star formation, protostars generate fast jets. The supersonic 
impact between the jet and the parent cloud generates shock fronts propagating through 
the high-density gas surrounding the protostar. The shocks heat and compress the ambient 
gas, and complex chemistry driven by endothermic chemical reactions, ice grain mantle 
sublimation and sputtering occur. Several molecules, such as H2O, NH3, CH3OH, H2CO, 
undergo abundance enhancements by orders of magnitude (e.g. Draine et al. 1983; Kaufman 
& Neufeld, 1996), as observed at mm-wavelengths in a number of outflows (Garay et al. 
1998; Bachiller & Perez Gutierrez 1997; Jimenez-Serra et al. 2005; J0rgensen et al. 2007). 

One of the best laboratories is the LI 157 region located at a distance between 250 
and 440 pc, where a Class protostar (L1157-mm) drives a chemically rich bipolar outflow 
(e.g, Gueth et al. 1996, Bachiller et al. 2001). The outflow is associated with several bow 
shocks seen in IR H2 (Nisini et al. 2010), in CO (Gueth et al. 1996) as well as in other 
molecules (Bachiller & Perez Gutierrez 1997) thought to be released from icy mantles such 
as H2CO, CH3OH, and typical tracers of high-speed shocks such as SiO, sputtered from dust 
grains (Caselli et al. 1997; Gusdorf et al. 2008; Jimenez-Serra et al. 2008). The brightest 
bow-shock, L1157-B1, has been extensively studied from the ground using IRAM-PdB and 
NRAO-VLA interferometers, reveahng a rich and clumpy structure (Tafalla & Bachiller 
1995; Benedettini et al. 2007; CodeUa et al. 2009). 

As part of the Herschel Key Program CHESS (Chemical Herschel Surveys of Star 
forming regions, Ceccarelli et al. 2010), L1157-B1 has been investigated with a spectral 
survey in the ~500-2000 GHz interval using the HIFI instrument to provide high-resolution 
(1 MHz) spectra. Preliminary results (CodeUa et al. 2010; Lefloch et al. 2010) have not 
only conflrmed the chemical richness of the LI 157- mm outflow, but have also revealed the 
presence of different molecular components at different excitation conditions coexisting in 
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the Bl structure. In addition, the CHESS data have allowed, for the first time, comparison 
of the line profiles of tracers of shocked material such as CH3OH, H2CO, with those due to 
the fundamental transitions of water and ammonia. These species show two different kinds 
of profiles, with CH3OH, H2CO, and NH3 emitting only at relatively low outfiow velocities 
(~15 km s~^ with respect to the ambient velocity), whereas H2O shows bright emission 
even at the highest velocities (up to 30 km s~^), as seen in SiO (e.g. Zhang et al. 1995) 
and CO (Lefloch et al. 2010). The intensity of methanol, formaldehyde, and ammonia with 
respect to H2O decreases towards higher outfiow velocities suggesting, in case of optically 
thin emission, a similar decrease in the abundance ratios. Note that while SiO and CO are 
not products of surface reactions, H2O, NH3, CH3OH and H2CO are all highly enhanced on 
the dust grains, and therefore the distinct behaviour of water is worth investigating. 

High NH3 abundances in the gas phase have been observed in many warm regions (e.g. 
hot cores and outfiows, Beuther et al. 2005; Sepulveda et al. 2011). NH3 is believed to be 
mainly a direct product of grain surface reactions (although at high temperatures ammonia 
can also form in gas phase, reaching abundances of ~ 10~^ - see Section 3) with nitrogen 
hydrogenating until saturation as it sticks to the dust during the cold phase (e.g. Brown et 
al. 1988; Viti et al. 2004a). Water is enhanced by the release of the icy mantles as well as 
by endothermic reactions occurring in the warm (> 220 K) shocked gas, which convert the 
bulk of atomic oxygen into water (e.g. Kaufman & Neufeld 1996, and references therein). 

In this Letter we investigate the differences between the water and ammonia profiles by 
using a new code which couples a gas-grain chemical model with a parametric shock model. 
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2. The model 

The gas-grain chemical code UCL_CHEM (Viti et al. 2004a) was coupled with the 
parametric shock model developed by Jimenez-Serra et al. (2008). UCL_CHEM is a time 
dependent gas-grain chemical model in two phases. During Phase I, gravitational collapse, 
gas-phase chemistry and sticking on to dust particles with subsequent processing occur. 
This phase simulates the formation of high density clumps and starts from a diffuse (~ 
100 cm~^) medium in atomic form (apart from a fraction of hydrogen in H2). The initial 
abundances adopted are solar but different initial values for the sulfur elemental abundance 
were also used. Atoms and molecules are depleted onto grain surfaces as in Rawlings et al. 
(1992, Eq. 2) and hydrogenate when possible. The depletion efficiency is determined by the 
fraction of the gas-phase material that is frozen onto the grains. This approach allows a 
derivation of the ice composition. Phase II computes the time dependent chemical evolution 
of the gas and dust once the clump has formed and stellar activity is present (in the form of 
a protostar and/or outflows). Full details of the code can be found in Viti et al. (2004a). For 
this work the UCL_CHEM was coupled with the parametric shock model of Jimenez-Serra 
et al. (2008) which calculates the physical structure of a plane-parallel steady-state C-shock 
that propagates with a velocity Vs through an unperturbed medium. This approximation 
was validated by comparing its results with those from magnetohydrodynamic (MHD) 
C-shock modeUing (Flower et al. 2003, Kaufman & Neufeld 1996). Although not strictly 
accurate, the general behaviour of the velocity of the ion and neutral fluids (i.e. Vi and 
Vn), of the density of the neutral particles (n„), and of the temperature of the ions and 
neutrals (Tj and r„), mimics the MHD physical structure of C-shocks (Figures 1 and 2 in 
Jimenez-Serra et al. 2008). 

In addition to the C-shock physical structure, the parametric model includes the 
sputtering of dust grains following the formalism used by Caselli et al. (1997), updated with 
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the results from May et al. (2000). In our chemical model, the icy mantles are sputtered 
once the dynamical age across the C-shock has reached the 'saturation time-scales' (or 
tsau Jimenez-Serra et al. 2008). These are the time-scales for which almost all molecular 
material within the mantles is released into the gas. In Jimenez-Serra et al. (2008), tsat 
were originally derived for SiO (see their Table 5). However, we adopt the same tsat values 
for II2O and NII3 because the mathematical expressions of the fraction of SiO, II2O and 
NH3 sputtered from the grain mantles are the same for all these molecular species except 
for their initial abundance contained in the mantles (see Appendix B in Jimenez-Serra et 
al. 2008). The adopted tgat for 112© and NH^ were finally scaled for the different initial H2 
densities explored in our study, by assuming [J that tgatocn'^. 

For this work. Phase I was used to form the pre-existing (to the shock) clump, Bl. 
We do not know how dense the pre-existing clump was; however modelling of the clumps 
along a similar outflow (CBS; Viti et al. 2004b; Benedettini et al. 2006) has shown that in 
order for the chemistry to be as rich as observed, some level of dumpiness must pre-exist 
the shock; hence we did not consider flnal densities for Phase I lower than 1000 cm^'^ (see 
Table 1). The shock model was then included in the Phase II of the chemical model. 

A total of 15 models were run (Table 1) where we varied: i) the final densities of Phase 
I (and hence the initial density of the clump before it became shocked); we have restricted 
these densities to the range between 10^ to 10^ cm~^, similar to those derived by Giisdorf 
et al. (2008) and Nisini et al. (2010) for the Bl clump. Higher pre-shock densities have 
not been considered in this study because the collisional dissociation of H2 would highly 
overcome the formation rate of H2, leading to the development of a J- type discontinuity 
within the physical structure of the shock for the shock velocities we are considering (Le 



^The cooling time-scales vary roughly as n- ^, where Ui is the density of the ion fiuid and 
is proportional to the H2 density of the gas - see e.g. Chieze et al. 1998; Lesaffre et al. 2004). 
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Bourlot et al. 2002). However, we note that, although most of the shocks along the 
LI 157- mm outflow are non-dissociative (or C-type; see Nisini et al. 2010), we cannot rule 
out the possibility that a compact J-type shock component is present in clump Bl, as 
suggested by the detection of [Sill], [Fell] or [SI] emission arising from this condensation 
(Neufeld et al. 2009). The flnal density and initial velocity also determine the maximum 
temperature attainable by the neutral gas within the shock (Draine et al. 1983 - see Figures 
8b and 9b for the assumed values of Tn,max shown in Table 1; Jimenez-Serra et al. 2008). 
Flower & Pineau des Forets (2003) predict higher temperatures for the same density and 
initial velocities at the beginning of the shock than Jimenez-Serra et al. (2008), although 
the values converge in the middle of the dissipation region. Since, however, recent MHD 
calculations from van Loo et al. (2009) seem to agree with the lower estimates, we adopt 
the calculations from Jimenez-Serra et al. (2008); ii) the initial elemental sulfur abundance: 
we use 1, 1/10 and 1/100 of the solar value; (iii) the degree of depletion, from 1 to 60% of 
gas depleted by the end of Phase I, an arbitrary range; (iv) the shock velocity from 35 to 
60 kms~^. 

Finally as a first attempt at reproducing the trends of the molecular line profiles, we 
input the fractional abundances, densities and temperatures as a function of distance and 
velocity produced by the chemical and shock model to the radiative transfer model SMMOL 
(Rawlings & Yates 2001), along with values for the dust temperature and radiation field 
intensity. The SMMOL code has accelerated A-iteration, a numerical scheme that solves the 
exact radiative transfer problem in multilevel non-local conditions, using an approximate 
lambda operator to derive 'pre-conditioned' statistical-equilibrium equations (e.g. Rybicki 
and Hummer 1991). We ran SMMOL for 0-H2O, and NH3. The coUisional data for water 
and ammonia were taken from Faure et al. (2007) and Danby al. (1988) respectively. New 
coUisional data for ammonia are available (Maret et al. 2009); however as they only cover 
up to temperatures of lOOK and para-H2, we have, for this work, used the older data. While 
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Table 1: Model parameters: Model number, pre-shock density, shock velocity, saturation 
time, maximum temperature of the neutral gas, degree of depletion at the end of Phase I, 
and initial sulfur abundance (where 1 stands for solar). 



Model 


n(H2) 


Vs 




T 

-■- n,max 


Depl. 


S 




cm~^ 


km s~^ 


yrs 


K 






1 


10^ 


40 


4.6 


4000 


1% 


1 


2 


10^ 


40 


4.6 


4000 


15% 


1 


3 


10^ 


40 


4.6 


4000 


30% 


1 


4 


10^ 


35 


4.6 


3200 


30% 


1 


5 


10^ 


35 


4.6 


3200 


1% 


1 


6 


10^ 


40 


455 


2200 


10% 


1 


7 


10^ 


40 


45.5 


2200 


15% 


1 


8 


10^ 


60 


38 


4000 


30% 


1 


9 


10^ 


60 


380 


4000 


30% 


1 


10 


10^ 


60 


380 


4000 


10% 


1 


11 


10^ 


40 


4.6 


4000 


30% 


1/10 


12 


10^ 


40 


4.6 


4000 


30% 


1/100 


13 


10^ 


40 


4.6 


4000 


60% 


1 


14 


10^ 


40 


45.5 


2200 


60% 


1 


15 


10^ 


35 


4.6 


3200 


60% 


1 
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water had been modelled by SMMOL before (Lerate et al. 2010) this is the first time that 
SMMOL was used to model the hne intensities and profiles of ammonia. This is the first 
attempt at using a non LTE non LVG radiative transfer code to model the high resolution 
H2O and NH3 hues observed by HIFI. 

3. Results 

The aim of this work is to understand qualitatively the origin and history of the 
molecular gas. Codella et al. (2010) find that the NH3/H2O intensity ratio decreases at high 
velocities implying a similar decrease in the abundance ratios. This is the first statement 
that we will test: Figure 1 shows the chemical evolution of selected species as a function of 
distance within the shock (z in units of cm) for four of our models (Model 1, 3, 4 and 7), 
chosen as examples of variations in depletion, pre-shock densities and shock velocities. We 
note that the value for the initial abundance of sulfur does not affect the water and ammonia 
behaviour at all and we won't discuss it further. An increase in depletion efficiency has the 
effect of increasing all the species that are formed or enhanced on the dust surfaces, but does 
not affect the abundance trends. Varying the pre-shock density and the shock velocity (and 
hence the neutral gas shock temperature) has a profound effect on the abundances trends. 
For Models 1 and 3 water behaves differently from the other three species in that it always 
remains abundant after the gas cools down, i.e it remains abundant at higher velocities 
within the post-shock gas (v^ > 15 km s~^). At first glance, this is certainly consistent with 
the assumption proposed by Codella et al. (2010) that the NH3/H2O decrease in intensity 
at high velocities corresponds to a decrease in their abundance ratio. In Models 4 and 7, on 
the other hand, ammonia remains high in abundance at high velocities. The key difference 
among these models is the maximum temperature reached by the gas due to the passage of 
the shock, and this, in turn, constrains the minimum pre-shock density for the clump (for 
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models with the same shock velocities). Indeed, if = 40 km s~^ only models with initial 
H2 densities greater or equal to 10^ cm~^ show temperatures of the neutral fluid of ~ 4000 
K in the postshock gas, reproducing the observed differences in the NH3/H2O abundance 
ratio. The fact that Model 7, where the initial density of the Bl clump is only 10^ cm~^, 
does not fit the observed behaviour is consistent with previous models of chemically rich 
clumps along outflows (Viti et al. 2004b), and supports the idea of an inhomogeneous 
medium with dense pre-existing clumps along L1157. Finally we also note that the observed 
decrease in the H2CO/H2O and CH3OH/H2O abundance ratios with increasing velocity 
within the shock (Codella et al. 2010), is better reproduced by Models 1 and 3 than by 
Models 4 and 7. 

For ammonia to be abundant in the gas phase alone, high temperatures are required; 
by running several UCL_CHEM Phase I test models (with constant densities of 10^ cm~^, 
with and without freeze out) we find that high abundances (at least 10~^) of pure gas phase 
ammonia are only reached if: (i) the temperature is higher than 800 K (as its main route 
of formation, H2 + NH2, is endothermic, with a barrier of ~ 1400K); (ii) no substantial 
freeze out occurs; (iii) the density of the gas is high (in our case 10^ cm~^) for at least 10^ 
years. These results are of course consistent with previous findings (see review by Ho & 
Townes 1983). Taking as an example Model 3, the imphcations for our LI 157 model are 
as follows: during Phase I of our model, when the temperature is lOK, the bulk of the 
ammonia is formed on the grains. During Phase II, the bulk of the ammonia is released 
from the grain mantles due to sputtering at tsat ~ 5 years (i.e at 2; ~ 5x10^^ cm); however 
an increase in abundance (from 4 to 6x10^^) does occur during the passage of the shock. 
The increase corresponds to when the temperature reaches ~1200K (corresponding to v„ ~ 
15 kms~^), but it is short-lived as NH3 starts dechning when the temperature reaches 3500K 
(corresponding to v„ ~ 25 kms~^). In terms of emission size of the clump this corresponds 
to less than 0.001 pc (0.5 -0.8 assuming a distance of 400 and 250 pc respectively). This 
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implies that observations of NH3 toward the Bl clump with single dishes (~20"-40") would 
then give lower NH3 abundances than those predicted by our model due to beam dilution 
(i.e. 5x10-5 X [l7(20"-40")2] - 10-^-3x10"^). We note that these diluted abundances 
are similar to those reported by Bachiller et al. (1993) toward the Bl clump. High-angular 
resolution observations of high-excitation hues of NH3 with the VLA are therefore required 
to establish the length of the region within the postshock gas where NH3 is expected to be 
enhanced. 

The degree and extent of the increase in the NH3 abundance due to gas phase high 
temperature chemistry will also depend on the degree of depletion adopted in Phase I; for 
example, if one compares Model 1 and Model 3 (where 1% and 30% of the gas phase is 
depleted by the end of Phase I, respectively) one notes that in Model 1 the relative increase 
is more substantial (at ~ 3x10^^ cm), due to the lower NH3 abundance released from the 
grains in Phase II. 

Ultimately, therefore, the difference between the water and ammonia behaviours 
between Models 1 and 3, and Models 4 and 7 is due to the temperature: in a model where 
the initial density is < 10^ cm~^, or the velocity is < 40 kms~^, the maximum temperature 
reached by the C-shock is ~ 3000 K. Under these conditions ammonia is maintained high 
as there is no effective destruction route while it can efficiently form via the reaction H2 
-I- NH2. For models where the temperature during the shock is ~ 4000 K, ammonia is 
instead efficiently destroyed by the reaction: H -|- NH3 NH2 + H2. The latter has a high 
barrier (~ 5000K) and hence it is inefficient for models where the maximum temperature of 
the shock is lower than 4000K. The main reaction that forms water, on the other hand, is 
very efficient as long as temperatures are over 220K while the main destruction route (via 
reaction with atomic hydrogen) has a very high barrier (~ 9000K). 

As the discriminating factor seems to be the maximum temperature attained by 
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the shock we note that if the shock velocity is higher (e.g. 60 kms~^) then low initial 
densities may still lead to maximum temperatures of ~ 4000K (Models 8-10). However 
we find that, regardless of the shock parameters, low initial densities imply that the gas 
abundances of water and ammonia (and other species) do not reach high values until the 
shock is at least half way through the dissipation region, leading to underestimated column 
densities. Moreover, if a shock velocity of 60 kms^^ were present toward Bl, we should 
expect to detect H2O at much higher terminal velocities than observed (~30 kms~^). These 
qualitative comparisons indicate that the best fit dynamical model for LI 157 is one where 
the shock velocity is ~ 40 kms~^ and n^^ ~ 10^ cm~^. 

Comparing our water abundances with those found by Lefloch et al. (2010) we 
are in good qualitative agreement: they find that the water is enhanced by two orders 
of magnitude between the low (from Ylsr='^-^ to -6 km s~^) and high-velocity (from 
Vi5R=-6 to -25 km s~^) components from 8xl0~^ to 8xl0~^ respectively; comparing these 
abundances with our lowest pre-shock water abundances of 8xl0~^ (corresponding to Ylsr 
up to -6 kms^^) and the high abundance reached during the shock of ~ 2xl0~^ (from -7 
to -30 kms~^), respectively; the latter is a factor of 3 higher than observed; we consider a 
match within a factor of three as a good agreement. 

Finally, we test whether the profile of the water line does indeed differ from that of 
ammonia. We have used the velocity, abundance, temperature and density profiles from 
Models 3 and 4 as input to our radiative transfer model SMMOL and we computed the hne 
profiles for water and ammonia. Figure 2 shows water and ammonia profiles for Model 3, 
overplotted with the ammonia profile for Model 4 (the water profile for Model 4 is very 
similar to that of Model 3). Clearly, while certainly not a match, the line profile for water is 
different from that of ammonia for Model 3, and in particular it extends to higher velocities. 
Moreover, in very good agreement with the observations, the water intensity is about a 
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factor of 7 higher than that of ammonia (the latter is multiphed by this factor in the figure 
- see caption) . On the contrary, the NH3 profile for Model 4 extends at higher velocities and 
its intensity is higher than observed (in the figure it is multiplied by a factor of 5) . This is 
the first time that modelling very high resolution profiles of water and ammonia in shocked 
regions starting from a pure ab initio approach has been attempted. Note that we do not 
aim at a quantitative match with the high resolution HIFI spectra: this would require a 
better knowledge of the geometry and structure of the Bl clump and most importantly a 
three dimensional radiative transfer model. Nevertheless this first attempt is very promising 
as we are in fact able to reproduce how the NH3/H2O intensity ratio changes with velocity. 
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Fig. 1. — Fractional abundances (with respect to molecular hydrogen), of selected species 
as a function of distance for Phase 11 for Models 1, 3, 4 and 7. The temperature profile for 
each model is also shown. 
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Fig. 2. — Theoretical line profiles of H20(lio-loi) (solid line), and NH3(lo-Oo) (dotted line) 
for Model 3. The profile of the NH3 hne has been multiphed by a factor of 7.5. The Vlsr is 
+2.6 kms~^. We also overplot the theoretical line profile of NH3(lo-0o) (dot-dash line) for 
Model 4, multiplied by a factor of 5. 
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